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Abstract 

We study a generalized Hadamard walk in one dimension with three inner states. The particle 
governed by the three-state quantum walk moves, in superposition, both to the left and to the 
right according to the inner state. In addition to these two degrees of freedom, it is allowed to 
stay at the same position. We calculate rigorously the wavefunction of the particle starting from 
the origin for any initial qubit state, and show the spatial distribution of probability of finding the 
particle. In contrast with the Hadamard walk with two inner states on a line, the probability of 
finding the particle at the origin does not converge to zero even after infinite time steps except 
special initial states. This implies that the particle is trapped near the origin after long time with 
high probability. 
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I. INTRODUCTION 

Quantum walks DDDDD are very simple quantum processes, however they are con- 
nected with a wide variety of field. Studies of the quantum walks are initially motivated by 
developing techniques for quantum algorithms. For example, the Grover's search algorithm 
[f|, which is one of the most famous quantum algorithms, is related to a discrete quantum 
walk ?l ■ Oka et al. |9j studied breakdown of an electric-field driven system by mapping 
the Landau-Zener transition dynamics to a quantum walk. 

The Hadamard walk plays a key role in studies of the quantum walk and it has been 
analyzed in detail. Thus the generalization of the Hadamard walk is one of fascinating 
challenges. The simplest classical random walker on a one-dimensional lattice moves to the 
left or to the right with probability 1/2. On the other hand, the quantum walker according 
to the Hadamard walk on a line moves both to the left and to the right. It is well known that 
the spatial distribution of probability of finding a particle governed by the Hadamard walk 
after long time is quite different from that of the classical random walk, see 



for examples. However there are some commonalities. In both walks the probability of 
finding probability at a fixed lattice site converges to zero after infinite long time. 

Let us consider a classical random walker in which it can stay at the same position with 
non-zero probability in a single time-evolution. In the limit of long time, dose the probability 
of finding a walker at a fixed position converges to a positive value ? If the probability of 
staying at the same position in a single step is very close to 1, the walker diffuses very 
slowly. However the probability of existence at a fixed position surely converges to zero if 
the jumping probability is not zero. The classical random walker who can remain the same 
position is essentially regarded as the same process with the random walk right or left with 
probability with 1/2 by scaling the time. Is the same conclusion valid for a quantum walk? 
The answer is "No" . We show that the profile of the Hadamard walk changes drastically by 
appending only one degree of freedom to the inner states. 

If the quantum particle in three-state quantum walk exists at only one site initially, the 
particle is trapped with high probability near the initial position. Similar localization has 
been already seen in other quantum walks. The first simulation showing the localization 
was presented by Mackay et al. Q| in studying the two-dimensional Grover walk. After 
that, more refined simulations were performed by Tregenna et al. and an exact proof 
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on the localization was given by Inui et al. 16]. The second is found in the four-state 

nn 

quantum walk [17, 18]. In this walk, a particle moves not only the nearest sites but also 
the second nearest sites according to the four inner states. The four-state quantum walk 
is also a generalized Hadamard walk, and it is similar with the three-state quantum walk. 
The significant difference between them is that the wavefunction of the four-state quantum 
walk dose not converge, but that of the three-state quantum walk converges in the limit 
as time tends to infinity. We focus on localized stationary distribution of the three-state 
quantum particle and calculate it rigorously. Moreover Konno proved that a weak limit 
distribution of a continuous-time rescaled two-state quantum walk has a similar form to that 
of the discrete-time one. In fact, both limit density functions for two-state quantum walks 
have two peaks at the two end points of the supports. As a corollary, it is easily shown that 
a weak limit distribution of the corresponding continuous-time three-state quantum walk 
has the same shape as that of the two-state walk. So the localization does not occur in the 
continuous-time case in contrast with the discrete-time case. In this situation the main aim 
of this paper is to show that the localization occurs for a discrete-time three-state quantum 
walk rigorously. 

The rest of the paper is organized as follows. After defining the three-state quantum 
walk, eigenvalues and eigenvectors of the time evolution operator are calculated, and the 
wave function is shown in Sec. II. A time-averaged probability of finding the particle is 
introduced in Sec. IE, and it is shown that the time-averaged probability of the three-state 
quantum walk converges to non-zero values. In Sec. IV, we prove the probability of finding 
the particle at a fixed itself converges to a non-zero value after infinite long time. 



II. DEFINITION OF THE THREE-STATE QUANTUM WALK 

The three-state quantum walk (3QW) considered here is a kind of the generalized 
Hadamard walk on a line. The particle ruled by 3QW is characterized in the Hilbert space 
which is defined by a direct product of a chirality state space \s) G {\L), |0), \R)} and a 
position space |n) G {. . . , | — 2), | — 1), |0), |1), |2), . . .}. The chirality states are transformed 
at each time step by the next unitary transformation: 

\L) = \{~\L) + 2|0) + 2\R)\ |0> = 1 -{2\L) - |0> + 2|i?», \R) = ±(2|L) + 2|0) - \R)). 
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Let ^(n, t) = [4>L(n, t), ipo(n, t), ipR^n, t)] T be the amplitude of the wave function of the 
particle corresponding to the chiralities "L", "0" and "R" at the position n £ Z and the 
time t £ {0, 1, 2, . . .}, where T denotes the transpose operator and Z is the set of integers. 
We assume that a particle exists initially at the origin. Then the initial quantum states are 
determined by [ip L (0, 0), ip (0, 0), ^(0, 0)] = [a, (3, 7], where a,/?, 7 £ C with |a| 2 + |/?| 2 + 
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2 = 1. Here C is the set of complex numbers. 
Before we define the time-evolution of the wavefunction, we introduce the following three 
operators: 



U, 
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If the matrix Ul is applied to the function ^(n, t), the only "L" -component is selected 
after carrying out the superimposition between ipL( n , t), ipo( n , t) and ipuin, t). Similarly the 
"0"-componet and "R"-componet are selected in Uo^(n,t) and Un^f(n,t). 
We now define the time evolution of the wavefunction by 

${n,t) = U L if>(n+l,t) + U if>(n,t) + U R il>(n-l,t). 

One finds clearly that the chiralities "L" and "R" are corresponding to the left and the right, 
and the chirality "0" is corresponding to the neutral state for the motion. 

Using the Fourier analysis, which is often used in the calculations of quantum walks, we 
obtain the wavefunction. The spatial Fourier transformation of ty(k,t) is defined by 

*(Jfe,t) = J2^{n,t)e~ ikn . 
The dynamics of wavefunction in the Fourier domain is given by 
















1 













e -lk 





1 

3 



-12 2 
2-12 
2 2-1 



*(M) 



(i) 



Thus the solution of Eq. (0) is formally given by ty(k,t) = U t ^(k, 0). Let e t9j - k and |$^) be 
the eigenvalues of U and the orthonormal eigenvector corresponding to e •*>* (j = 1,2,3). 
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Since the matrix U is a unitary matrix, it is diagonalizable. Therefore the wavefunction 
ty(k,t) is expressed by 

where \f(/c, 0) = [a,/?,7] T G C 3 with |a| 2 + |/?| 2 + |7| 2 = 1. The eigenvalues are given by 



cos #i 



0, j = l, 
k , 3 = 1, 
-0 k , 3 = 3, 

-i(2 + cosA;), 



sm 



k = -V (5 + cos k)(l - cos/c), 



(2) 



for /c G [— 7r, 7r). The eigenvectors of U which are orthonormal basis are obtained after some 
algebra: 



c k (0 



j,k) 



where 



c k (0) 
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(3) 



1 + cos(6» - k) 1 + cos 9 1 + cos(6 l + k) 
In the above calculation, it is most important to emphasize that there is the eigenvalue 
"1" independently on the value of k. The eigenvalues of the Hadamard walk are given by 
e t6k and e 1 ^ 9 ^ with the arguments satisfying sin^ = smk/y/2. Therefore the eigenvalues 
does not take the value "1" except the special case k = 0. We have shown the existence 
of strongly degenerated eigenvalue such as "1" in Eq. (J2J) is necessary condition in the 
quantum walks showing the localization Q]. The significant difference between 3QW and 
the previous quantum walks showing the localization such as the Grover walk and the four- 
state quantum walk is the number of the degenerate eigenvalues of the time evolution matrix. 
There are two degenerate eigenvalues with values "1" and independently of the value 
of k in both the Grover walk and the four-state quantum walk. In contrast there is only 
one degenerate eigenvalue independently of the value k for 3QW. This distinctive property 
causes the particular time evolution in 3QW. 



Let us number each chirality "L","0", and "R" using I = 1, 2, and 3, respectively. Then 
the wavefunction in real space is obtained by the inverse Fourier transform: for a,/3, 7 G C 
with \a\ 2 + \(3\ 2 + |7| 2 = 1, 

tf(n,t;a,/3,7) = — [ ^(k,t)e ikn dk 
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where 



0=1 

[*(n, *; 1; a, /?, 7), tt(n, t; 2; a, /?, 7), t; 3; a, /3, 7 )] T 
3 

^ [%(n, t; 1; a, /3, 7), ^(n, t; 2; a, /?, 7), tf^n, t; 3; a, /3, 7 )] T , (4) 
i=i 



^■(n,t;Z;a,/?, 7 ) = ^ / c k (O jtk )<p k (0 j>k , l)e l ^ t+kn ^dk, (Z = 1,2,3) (5) 



with 



^(0, = Cz,fc(0)[«Ci,fc(0) + /^(G) + 7C 3 , fc (^)], (J = 1, 2, 3) 

Ci,^) = (l + e l(e - fc) r\ C2, fc W = (l + e t T\ CM(e) = (l + e'(* + * ) )- 1 , (6) 

where z is conjugate of z G C. Sometimes we omit the initial qubit state [a, (3, 7] such as 
tyj(n, t; I) = ^j{n, t; I; a, /3, 7). The probability of finding the particle at the position n and 
time t with the chirality I is given by P(n, t; I) = \^(n, t; l)\ 2 . Thus the probability of finding 
the particle at the position n and time t is P(n, t) = Ya=i P(. n i t\ 0- 



III. TIME-AVERAGED PROBABILITY 



We focus our attention on the spatial distribution of the probability of finding the particle 
after long time. Equation (£Q) is rather complicated, but the value of lim^oo P(n, t) can 
be exactly calculated in following sections. We start showing a numerical result for the 
probability of finding the particle at the origin before carrying analytical calculation. Figure 
1 shows the time dependence of P(0,t) with the initial state a = i/y/2,/3 = 0,7 = l/y/2. 
The probability decreases quickly near t — and fluctuates near 0.2. The inserted figure 
in Fig. 1 is plotted as the inversed time, and the amplitude of fluctuation decreases as the 
time. The horizontal dashed line in Fig. 1 is the value at the origin of the time-averaged 
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probability defined by 

a, /3, 7) = lim I Jim ^y]y] Pn{0, t; I; a, (3, 7) 



/ 3 T-l 

lim lim - > > P^fO, t; I; 

\ 1=1 t=0 



where Pv(n, t; Z; a, (3, 7) is the probability of finding a particle with the chirality I at the 
position n and time t on a cyclic lattice containing N sites. The time-averaged probability 
introduced here has already used to study the quantum walk. In the case of the Hadamard 
walk on a cycle containing odd sites, the time-averaged probability takes a value 1/N in- 
dependently of the initial states jj, ■ Therefore the time-averaged probability converges 
to zero in the limit iV — > 00. On the other hand the time-averaged probability of quantum 
walks which exhibits the localization converges to a non-zero value. For this reason we firstly 
calculate time-averaged probability of 3QW, and we will show in the next section that the 
probability P(n,t) itself converges as well. 

We calculate the time-averaged probability at the origin in 3QW on a cycle with iV sites. 
We assume that the number N is odd. The argument of eigenvalues of 3QW with a finite 
N is Oj^nm/N f° r m £ [ — (N — l)/2, (N — l)/2]. Since the eigenvalues corresponding to m is 
the same with the eigenvalues corresponding to —m, the wavefunction is formally expressed 
by 

3 (JV-l)/2 
j=l m=0 

where ^jvj.m is Qj,2mw/N- We note here that the coefficients Cj !m< i(N) depend on the initial 
state [a, j3, 7], but we omit to describe it. Then the probability Pv(0, £; a, j3, 7) at the origin 
is given by 

3 (N-l)/2 
h,h,ji,h=l mi,m 2 =0 

The coefficients Cj tTn j(N) are determined from the product of eigenvectors. Although lengthy 
calculations are required to express the coefficients Cj, m ,i(N) explicitly, it is shown below 
that coefficients Cj iTrii i(N) except j — 1 dose not contribute Pjv(t). Noting that the following 
equation 



T-l 

lim ±Ye m 



t=o 




we have 



P N (0;a,(3,j) = 



1=1 



(N-l)/2 

E C ^m,l( N ) 



m=0 



i=2 



3 (N-l)/2 

E E i^wWi 1 



(7) 



j'=2 m=l 

The difference between the first term and the other terms is caused by the difference of the 
degree of degenerate eigenvalues. 

Let (pj t k{N) be the eigenvectors corresponding to the eigenvalue e l9j - kt of the time-evolution 
matrix of the 3QW with a matrix size 3N x 3N. They are easily obtained from |$^,) by 

1 



7 3iV 



where ui = e 2m ^ N . Since the coefficients Cj im> i(N) are proportion to the product of eigenvec- 
tors, the orders with respect to iV of the first term and the second term in Eq. are 0(1) 
and 0{N~ X ) ) respectively. Thus we can neglect the second term in the limit of N — > oo. 
Using the eigenvectors in Eq. Q, we have 

Pnfra,?^) = (5-2v / 6)(l + |a + /5| 2 + |/5 + 7 | 2 -2|/3| 2 ). 

The time-averaged probability takes the maximum value 2(5 — 2y / 6) at (3 — 0, which is the 
value indicated by the horizontal dashed line in Fig. 1. The component of -Poo(0; a, (3, 7) 
corresponding to I = 1, 2, 3 are respectively given by 

\y/Ea - 2( v / 6 - 3)/9 + (12 - 5V6)7| 2 



P^O; 1; a, /3, 7) 
£«,((); 2; a, 0, 7) 
3; a, 0, 7) 



36 



(y / 6-3) 2 |a + /j + 7 | 2 
9 

IV67 - 2(VQ - 3)(3 + (12 - 5^) 



a 



36 



(8) 



We stress here that the time-averaged probability is not always positive. If a = l/y/E, 
(3 = —2/^/6 and 7 = 1/y/E, then the time-averaged probability becomes zero, that is, 
PootO; 1/V6, -2/V6, 1/V6) = 0. 



IV. STATIONARY DISTRIBUTION OF THE PARTICLE 



We showed that the time- averaged probability of 3QW converges to non-zero values ex- 
cept special initial states. This result, however, does not mean that a particle is observed with 
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almost the same probability at the origin after long time evolution. Indeed the probabilities 
of finding a particle in Grover walk and the four-state quantum walk, whose time-averaged 
probability converge to non-zero values, do not converge. Because the time-averaged proba- 
bilities in both cases depend on the parity of the time. Figure 1 suggests that the probability 
P(n,t) = P(n,t;a, ft,j) itself converges in a limit of t — > oo. In this section we calculate 
the limit P*(n) = limbec P(n, t) rigorously and consider the dependence of P*(n) on the 
position n. 

The wavefuntion given by Eq. (J3J) is infinite superimposition of the wavefunction 
e i(6j,kt+ kn ) . If the argument Oj^ given in Eq. (J2J) for j = 2,3 is not absolute zero, 
then the e t0 ^ kt oscillates with high frequency in /c-space for large t. On the other hand, 
Ck(9j,k)(Pk(@j,kJ) in Eq. (jSJ) changes smoothly with respect to k. Therefore we expect that 
the integration for j = 2, 3 in Eq. (jSJ) becomes small as the time increases due to cancellation 
and converges to zero in the limit of t — > oo. That is, 



for / = 1,2,3 and a,/?, 7 G C with \a\ 2 + \ft\ 2 + I7I 2 = 1. This conjecture can be rigorously 
proved using the Riemann-Lebesgue lemma (see Appendix A). Consequently the probability 
P*(ti) = P*(n; a,/?, 7) is determined from the eigienvectors corresponding to the eigenvalue 
"1", that is, 9 l k = (see Eq. ((2))), and l-th component of P*(n; I) = P*(n; I; a, j3, 7) is given 



Note that \I/i(n, t; I; a, ft, 7) does not depend on time t, since 61 k = 0. By transforming 
integration in the left-hand side in Eq. (fTU|) into complex integral, we have 



3 




(9) 



by 



P m (n;l;a,ft,i) = |*i(n, t; Z; a, (3, 7)| 2 . 



(10) 



P*(n;3;a,/3,7) 



P*(n; l;a,ft,j) 



P*(n;2;a,/3,7) 



\2al(n) + ftJ+(n) + 2 1 K + (n)\ 2 , 

ft 2 
aJ-{n) + -^L{n) + ^J+{n) , 

\2aK_(n) + ftj_(n) + 2 7 /(n)| 2 , 



(11) 



where c 



5 + 2^(6 (-1,0)) and 
2c l n l +1 



I(n) = -5 — r , L(n) = /(n - 1) + 2I(n) + J(n + 1) 
c z — 1 

J+(n) = J(n) + /(n + 1), J_(n) = /(n - 1) + /(n), 



K + (n) = /(n + l), X_(n) =J(n-l), 



(12) 
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for any n G Z. We should remark that it is confirmed Px>(0; I; a, /3, 7) = P*(0; Z; a, /?, 7) for 
/ = 1, 2, 3 by using Eqs. ©, (fTT|) and (fT2j) . 

Here we give an example. From Eqs. and (|12j). we obtain 

P,(0; i/y/2, 0, l/y/2) = + ^ 2 + 5) = 10 - 4^6 = 0.202 . . . , (13) 

1 — c z Y 



P.dnM/^O,!/^) = 2(5^ + 2^lg + 2c + 5) c2 ,„, 



(14) 



for any n > 1. Furthermore, 

< ^P*(n;z/V2,0,l/V2) = 1/^6 = 0.408 .. . < 1. (15) 

That is, P*(n; i/y/2, 0, l/y/2) is not a probability measure. The above value depends on the 
initial qubit state, for example, 

P* {n; 1/V3, l/y/3, l/y/3) = 3 - y/6 = 0.550 . . . , 

n&L 

p *{n\ 1/V3, -1/V3, l/y/3) = (3 - Vq)/9 = 0.061 .... 

We should remark that in the case of the classical symmetric random walk starting from 
the origin, it is known that P*(n) = for any n G Z, therefore we have Yl n ezP*( n ) = 0- 
The same conclusion can be obtained for the discrete-time and continuous-time two-state 



quantum walks 
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Figure 3 shows the probability P* (n; i/y/2, 0,1/ \/2) in log scale. The probability decreases 
exponentially for large |n| and its asymptotic behavior is express by P*(n; i/y/2, 0, l/y/2) oc 
the limit of n — > ±cxo as Eq. (|14j) indicates. 

We here mention the time-dependence of the P(n,t). The particle is observed near the 
origin with high probability. This, however, dose not imply that the particle can not escape 
form the region near the origin. Figure 4 shows the change of P(n, t) in a space-time. The 
particle is observed with high-probability at the dark region. One clearly finds three dark 
regions. The first is a region near a center line connecting with the origin, and we can 
confirm the localization near the origin. The second and third regions are boundaries of the 
triangle. These regions show trajectories of two peaks moving outside in the space-time, 
which are also seen in the Hadamard walk. As a result we conclude that the particle in 
3QW splits three parts in superposition. 
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V. CONCLUSIONS AND DISCUSSIONS 



The unique properties which are not observed in other quantum walks were found in 
three-state quantum walk. The particle which exists at the origin splits three pieces in 
superposition. Two of three leave for infinite points and the remainder stays at the origin. 
In contrast with the ordinary Hadamard walk, the probability of finding the particle at the 
origin does not vanish for large time, and its maximum probability is 10 - 4^ = 0.202 .... 

A simple reason why the 3QW is different form other quantum walks is the difference 
in the degree of degenerate eigenvalues. The necessary condition of the localization is the 
existence of the degenerate eigenvalues. And furthermore, the each degree of degeneration 
must be proportion to the dimension of the Hilbert space. In addition, if the degenerate 
eigienvalue is 1 only, then the probability of finding the particle can converges in the limit of 
t — > oo. The Grover walk and the four-state quantum walk exhibit the localization, but the 
probability of finding the particle oscillates. Because there are degenerate eigenvalues with 
values "1" and "-1". On the other hand, the degenerated eigenvalue in 3QW unrelated to 
the wave number is only "1" , therefore the probability converges. Although no experiment 
exists about the quantum walks yet, the stationary properties of 3QW may be advantage to 
comparison with theoretical results. 

Finally we discuss a relation between the limit distribution for the original 3QW X t as 
time t — > oo and that of the rescaled X t /t in the same limit. When we consider the 3QW 
starting from a mixture of three pure states [1, 0, 0] T , [0, 1, 0] T , and [0, 0, 1] T with probability 
1/3 respectively, we can obtain a weak limit probability distribution f(x) for the rescaled 
3QW Xt/t as t — > oo in the following: 

for x G [—1,1], where 5 (x) denotes the pointmass at the origin and L at b)(x) — 1, if x 6 



12 



(a, 6), = 0, otherwise. The above derivation is due to the method by Grimmett et al 
We should note that the first term in the right-hand side of Eq. (fT7)|) corresponds to the 
localization for the original 3QW. In fact, from Eq. ()12j) . we have 

~ l p M 1, 0, 0) + P„(n; 0, 1, 0) + P*(n; 0, 0, 1)] = 1 

The last value 1/3 is nothing but the coefficient 1/3 of 5q(x). Moreover, the second term in 
the right-hand side of Eq. (JTfij) has a similar form a weak limit density function for the same 
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rescaled Hadamard walk with two inner states starting from a uniform random mixture of 
two pure states [1, 0] T and [0, 1] T : 

f / \ = / (-l/v / 2,l/v / 2)( X ) 

1 ] 7r(l-x 2 )v / l-2x 2 ' 

for x G [—1,1], see jlO, 11, 3|. This case does not have a delta measure term corresponding 
to a localization. More detailed study on this line will appear in our forthcoming paper. 



APPENDIX A: PROOF OF EQ. © 

Here we give an outline of the proof of Eq. ©. A direct computation gives 
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*i 


(n,*; 1; 


a,/3,7) 
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*i 


(n,t;2; 




= M 


/? 


i=2 


*i 


(n,*;3; 


a,/3,7) 
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where M = (mi,)i<j 7 o with 



mn — 3J ni t + — { J n -\,t + Jn+l,t + 

m 33 — 3J ni t + - {J n -l,t + Jn+l,t — (Kn-l,t 

m 12 — — \ Jn,t + Jn+l,t + (-^n,t — ^n+l,t)} j 

™32 = - {Jn,t + Jn-l,t + {K n ,t - K n-l,t)} , 

mi3 — ~^Jn+l,t, ^31 = — 2J n _i >t , 

m 21 = — {Jn,t + Jn-l,t + (-^n-l,t — ^n,t)} ? 

^23 — — {Jn,t + Jn+l,t + (-^n+l,t — ^n,t)} ? 



™22 = 4J n 



and 



•^n.t — - — 



" cos(fcn) 
27T J_ w 5 + cos 



cos(9 k t)dk, 



K n+ i,t)} 



K n , t = — 



cos(kn) 



27r J_ n a/(5 + cos k)(l - cos A;) 
From the Riemann-Lebesgue lemma, we can show that 



sm(9kt)dk. 



lim J n t = 0, 

t— >oo 



lim (K n t - K n+ i )t ) 

t— >oo 



0. 
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for any nsZ. Therefore we have the desired conclusion. 
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FIG. 1: Time dependence of the probability of finding a particle at the origin starting from an 
initial state a = i/\/2,(3 = 0, 7 = l/\/2. The same data are plotted as a function of 1/t in the 
inserted figure. The horizontal dashed line cross at (0,2(5 — 2\/6)). The amplitude of fluctuation 
of probability decreases near zero. 
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FIG. 2: Distribution of probability of finding a particle P*(n) at the position n in 3QW with 
the initial state a = i/V2,/3 = and 7 = 1/V2 on a semi- log scale. The probability decreases 
exponentially for large \n\. 
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FIG. 3: Density plot of probability of finding a particle P(n, t) in the space-time. The initial state 
is a = i/y/2,/3 = 0,7 = 1/V2. The particle is observed on darker regions with high probability. 
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